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Coherent large-scale circulations of turbulent thermal convection in air have been studied ex- 
perimentally in a rectangular box heated from below and cooled from above using Particle Image 
Velocimetry. The hysteresis phenomenon in turbulent convection was found by varying the temper- 
ature difference between the bottom and the top walls of the chamber (the Rayleigh number was 
changed within the range of lO'' - 10*). The hysteresis loop comprises the one-cell and two-cells 
flow patterns while the aspect ratio is kept constant {A = 2 — 2.23). We found that the change of 
the sign of the degree of the anisotropy of turbulence was accompanied by the change of the flow 
pattern. The developed theory of coherent structures in turbulent convection (Elperin et al. 2002; 
2005) is in agreement with the experimental observations. The observed coherent structures are 
superimposed on a small-scale turbulent convection. The redistribution of the turbulent heat flux 
plays a crucial role in the formation of coherent large-scale circulations in turbulent convection. 



I. INTRODUCTION 

Coherent structures in turbulent convection are ob- 
served in the atmospheric convective boundary layers 
(see, e.g., Etling and Brown 1993; Atkinson and Wu 
Zhang 1996; Briimmer 1999) and in laboratory exper- 
iments (see, e.g., Krishnamurti and Howard 1981; Sig- 
gia 1994; Kadanoff 2001). In the atmospheric shear-free 
convection, the coherent structures represent large-scale 
three-dimensional long-lived Benard-type cells (cloud 
cells) composed of narrow uprising plumes and wide 
downdraughts. They usually embrace the entire convec- 
tive boundary layer (of the order of 1-3 km in height) 
and include pronounced convergence flow patterns close 
to the surface. In the sheared convective flows, the coher- 
ent structures represent convective boundary layer scale 
rolls (cloud streets) stretched along the mean wind (see, 
e.g., Etling and Brown 1993; Atkinson and Wu Zhang 
1996; Briimmer 1999). 

Buoyancy-driven structures, such as plumes, jets, and 
large-scale circulation patterns were observed in numer- 
ous laboratory experiments. The large-scale circulation 
caused by convection in the Rayleigh-Benard apparatus 
is often called the "mean wind" (see, e.g., Krishnamurti 
and Howard 1981; Zocchi et al. 1990; Ciliberto et al. 
1996; Niemela et al. 2001; Sreenivasan et al. 2002; 
Niemela and Sreenivasan 2003; Burr et al. 2003; Shang 
et al. 2004; Xi et al. 2004; Tsuji et al. 2005; Brown et 
al. 2005, and references therein). There are several open 
questions concerning these flows, e.g., how do they arise, 
and what are their characteristics and dynamics. 

The life-times of coherent structures are very long com- 
pared to the largest turbulent time-scales. Their spec- 
tral properties differ from those of small-scale turbulence. 
They are characterized by narrow spectra and do not ex- 
hibit the direct energy-cascade behavior (from larger to 
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smaller scales). As a result the turbulence and the co- 
herent structures interact in practically the same way as 
the turbulence and the mean flow. These structures show 
more similarity in their behavior with regular flows than 
with turbulence. They can be identified as the motions, 
whose spatial and temporal scales are much larger than 
the characteristic turbulent scales. 

In spite of a number of theoretical and numerical stud- 
ies (see, e.g., Busse and Whitehead 1971; 1974; Busse 
1983; Lenschow and Stephens 1980; Hunt 1984; Hunt 
et al. 1988; Schmidt and Schumann 1989; Zilitinkevich 
1991; Wilhams and Hacker 1992; 1993; Zihtinkevich et 
al. 1998; Young et al. 2002; Parodi et al. 2004, and ref- 
erences therein) , the nature of large-scale coherent struc- 
tures in turbulent convection is a subject of discussions. 
Hartlep et al. (2003) noted that there are two points 
of view on the origin of large-scale circulation in turbu- 
lent convection. "According to one point of view, the 
rolls which develop at low Rayleigh numbers, Ra, near 
the onset of convection continually increase their size as 
Ra is increased and continue to exist in an average sense 
at even the highest Rayleigh numbers reached in the ex- 
periments (see, e.g., Fitzjarrald 1976). Another hypoth- 
esis holds that the large-scale circulation is a genuine 
high Rayleigh number effect (see, e.g., Krishnamurti and 
Howard 1981)." 

The new mechanism of formation of the large-scale co- 
herent structures in turbulent convection was proposed 
recently by Elperin et al. (2002; 2005). It was suggested 
that the redistribution of the turbulent heat flux plays 
a crucial role in the formation of the large-scale circu- 
lations in turbulent convection. In particular, two com- 
petitive effects, namely redistribution of the vertical tur- 
bulent heat flux due to convergence or divergence of the 
horizontal mean flows, and production of the horizontal 
component of the turbulent heat flux due to the interac- 
tion of the mean vorticity with the vertical component of 
the turbulent heat flux, cause the large-scale instability 
and formation of the large-scale coherent structures in 
turbulent convection (Elperin et al. 2002; 2005). 
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The main goal of this paper is to describe the experi- 
mental study of large-scale circulations of turbulent ther- 
mal convection in air flow (the aspect ratios A = 2—2.23). 
In order to study large-scale circulations we used Parti- 
cle Image Velocimetry to determine the turbulent and 
mean velocity fields, and a specially designed tempera- 
ture probe with twelve sensitive thermocouples was em- 
ployed to measure the temperature field. We found the 
hysteresis phenomenon in turbulent convection by vary- 
ing the temperature difference between the bottom and 
the top walls of the chamber (the Rayleigh number was 
changed within the range of Ra = 10^ — 10^). The hys- 
teresis loop comprises the one-cell and two-cell flow pat- 
terns. We found that the change of the sign of the de- 
gree of the anisotropy of turbulence was accompanied by 
the change of the pattern of the mean wind. The devel- 
oped theory of coherent structures in turbulent convec- 
tion (Elperin et al. 2002; 2005) is in a good agreement 
with the experimental observations. 

The hysteresis phenomenon in laminar convection was 
found by Busse (1967) who defined it as follows: "The 
fact that the convection at a certain Rayleigh number 
depends on the way in which the Rayleigh number has 
been reached is called the hysteresis effect in laminar con- 
vection". In the laminar convection the hexagon flow 
structures transform into rolls structures by increasing 
the Rayleigh number. On the other hand, decreasing 
the Rayleigh number causes the transition from the rolls 
structure to the hexagons. Similar phenomena were ob- 
served in numerical simulations and laboratory experi- 
ments (see, e.g., Busse 1978; Braunsfurth et al. 1996; 
Gelfgat et al. 1999, and references therein). 

In the experiments by Willis et al. (1972) performed 
in air, water and a silicon oil for the Rayleigh numbers 
Ra = (0.2 — 3) X 10^ and large aspect ratios (for air 
^ = 31.5), it was demonstrated that the average dimen- 
sionless roll diameter, A, increases as Ra is increased. 
This observation was less pronounced for large Rayleigh 
number flows (for water and silicon oil). The hysteresis 
phenomenon was found in dependence A(Ra) for large 
Prandtl number flows (for water and silicon oil), but this 
phenomenon was not observed in the air flow (Willis et 
al. 1972). For considerably larger values of the Rayleigh 
numbers the convection patterns become very compli- 
cated (Willis et al. 1972) and turbulent convection arises. 
In the experiments by Fitzjarrald (1976) in air flow it 
was shown that in convection with large aspect ratios 
the predominant horizontal scale of the regular struc- 
tures increased from 4/i at Ra = 10^ to 6h at Ra = 10^ 
(where h is the vertical distance between plates). In the 
range Ra 10^ — 10^ this scale did not change (Fitzjar- 
rald 1976). Notably, the hysteresis phenomenon was not 
observed by Fitzjarrald (1976). 

In the direct numerical simulations of convection by 
Hartlep et al. (2003) for the Rayleigh numbers Ra = 
10^ — 10^, the Prandtl number Pr = 0.7 and large as- 
pect ratio {A = 10) it was found that "the typical size 
of the large-scale structures does not always vary mono- 



tonically as a function of Ra, but broadly increases with 
increasing Ra. It cannot be decided from these simula- 
tions whether the large-scale structures will eventually 
disappear at yet higher Rayleigh numbers." Notably, the 
hysteresis phenomenon was not observed by Hartlep et 
al. (2003). 

The paper is organized as follows. Section 2 discusses 
the physics of formation of the large-scale coherent struc- 
tures, and Section 3 describes the experimental set-up for 
a laboratory study of this effect. The experimental re- 
sults and their detailed analysis are presented in Section 
4. Finally, conclusions are drawn in Section 5. 



II. TURBULENT HEAT FLUX AND 
COHERENT STRUCTURES 

Traditional theoretical models of the boundary-layer 
turbulence, such as the Kolmogorov-type local closures, 
imply the following assumptions. Fluid flows are decom- 
posed into two components of principally different na- 
ture, organized mean flow and turbulent flow. Turbu- 
lent fluxes are proportional to the local mean gradients, 
whereas the proportionality coefficients (eddy viscosity, 
turbulent diffusivity) are uniquely determined by local 
turbulent parameters. For example, a widely used tra- 
ditional approximation for the turbulent heat flux reads 
F = (su) = -ktVS" (sec, e.g., Monin & Yaglom 1975), 
where kt is the turbulent thermal conductivity, S is the 
mean entropy, u and s are fluctuations of the velocity 
and entropy, respectively. 

The traditional form for the turbulent heat flux F does 
not include the contribution from anisotropic velocity 
fluctuations. Actually the mean velocity gradients can 
directly affect the turbulent heat flux. The reason is that 
additional essentially anisotropic velocity fluctuations are 
generated by tangling the mean-velocity gradients with 
the Kolmogorov-type turbulence due to the influence of 
the inertial forces during the life time of large turbulent 
eddies. Therefore, the Kolmogorov turbulence supplies 
energy to the anisotropic (tangling) turbulence. In its 
turn the anisotropic turbulence causes formation of co- 
herent structures due to the excitation of a large-scale 
instability (Elperin et al. 2002; 2005). Anisotropic veloc- 
ity fluctuations were studied for the first time by Lumley 
(1967). He had shown that the velocity field in the pres- 
ence of mean shear is strongly anisotropic and is charac- 
terized by a steeper spectrum than the Kolmogorov tur- 
bulence (see also Wyngaard and Cote 1972; Saddoughi 
et al. 1994; Ishihara et al. 2002; Yoshida et al. 2003, 
among others). 

In order to parameterize the anisotropic turbulence, a 
spectral closure model was developed by Elperin et al. 
(2002; 2005). The derivation of the anisotropic turbu- 
lence model includes the following steps: applying the 
spectral closure, solving the equations for the second 
moments in the k space, and returning to the physical 
space to obtain formulas for the Reynolds stresses and 
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the turbulent heat flux. The derivation are based on the 
Navier-Stokes equation and the entropy evolution equa- 
tion formulated in the Boussinesq approximation. This 
derivation yields the following expression for the turbu- 
lent flux of entropy F = (su): 



aFZ divU,, 



(WxF* 



(1) 



(for details, see Elperin et al. 2002; 2005), where r is 
the correlation time of the Kolmogorov turbulence cor- 
responding to the maximum scale of turbulent motions, 
W — VxU is the mean vorticity, U = XJh + is 
the mean velocity with the horizontal U/j and vertical 
components, a is the degree of thermal anisotropy, 



F* = 



, V j S is the background turbulent heat flux. 



(2) 



is a generalized anisotropic turbulent heat conductivity 
tensor, 7 is the ratio of specific heats and e is the vertical 
unit vector. The equation for the tensor Kij was derived 
in Appendix A in Elperin et al. (2002) using the budget 
equations for the turbulent kinetic energy, fluctuations of 
the entropy and the turbulent heat flux. The anisotropic 
part of the tensor Kij (the second term in the square 
brackets) is caused by a modification of the turbulent 
heat flux by the buoyancy effects. The parameter a is 
given by 



l + C/3 ' 
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where Ih and Iz are the horizontal and vertical scales 
in which the background turbulent heat flux F*{r) = 
(s(x) Uz(x-l-r)) tends to zero. The parameter ^ describes 
the degree of thermal anisotropy. In particular, when 
Ih = Iz the parameter ^ = and a — 1. For Ih <^ Iz the 
parameter ^ = —1 and a — —9/2. The maximum value 
^max of the parameter ^ is given by ^,nax = 2/3 for a 3. 
The upper limit for the parameter ^ arises because the 
function F*{r) has a global maximum at r = 0. Thus, for 
a < 1 the thermal structures have the form of columns 
or thermal jets (Ih < h), and for a > 1 there exist the 
"pancake" thermal structures {Ih > Iz) in the background 
turbulent convection (i.e., a turbulent convection with a 
zero gradients of the mean velocity). 

The terms in the square brackets in the right hand 
side of Eq. are caused by the anisotropic turbulence 
and depend on the "mean" (including coherent) velocity 
gradients. These terms lead to the excitation of large- 
scale instability and formation of coherent structures. In 
Eq. (^) the terms with zero divergence are omitted, be- 
cause only divF contributes to the mean-field dynamics. 
Neglecting the anisotropic turbulence term this equation 
reduces to the traditional equation. The physical mean- 
ing of Eq. is the following. The second term in Eq. 
describes the redistribution of the vertical background 





FIG. 1: Effect of a nonzero div Uh which causes a redistribu- 
tion of the vertical turbulent flux of the entropy and results 
in a formation of a large-scale circulation of the velocity field: 
(a) fluid flow with div \Jh < produces regions with enhanced 
vertical fluxes of entropy and vertical fluid flow in these re- 
gions (b). 



turbulent heat flux by convergent (or divergent) horizon- 
tal mean flows. This redistribution of the vertical tur- 
bulent heat flux occurs during the life-time of turbulent 
eddies. The last term in Eq. ^ determines the forma- 
tion of the horizontal turbulent heat flux due to "rota- 
tion" of the vertical background turbulent heat flux by 
the perturbations of the horizontal mean vorticity. The 
emergence of the horizontal turbulent heat flux is caused 
by the action of the local inertial forces in inhomogeneous 
mean flows. 

In the shear-free regime, the large-scale instability is 
related to the first term in square brackets in Eq. ^ for 
the turbulent flux of entropy. When dUz/dz > 0, per- 
turbations of the vertical velocity Uz cause negative di- 
vergence of the horizontal velocity, div \Jh < (provided 
that divU = 0). This strengthens the local vertical tur- 
bulent flux of entropy and causes increase of the local 
"mean" entropy and buoyancy (see Fig. 1). The latter 
enhances the local "mean" vertical velocity Uz- Thus 
a large-scale instability is excited. Similar reasoning is 
valid when dUz/dz < 0, whereas divU/i > 0. Then 
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a negative perturbation of the vertical flux of entropy 
leads to a decrease of the " mean" entropy and buoyancy, 
that enhances the downward flow and once again excites 
the instability. Therefore, nonzero div \Jh causes redis- 
tribution of the vertical turbulent flux of entropy and 
formation of regions with large values of this flux. These 
regions (where divUji < 0) alternate with the low-flux 
regions (where divU^ > 0). By this means large-scale 
coherent structures are formed. 

The role of the second term in square brackets in 
Eq. ^ is to decrease the growth rate of the large-scale 
instability for a > —3/2. Indeed, the interaction of 
the "mean" vorticity with the vertical flux of entropy 
produces the horizontal turbulent heat flux. The latter 
decreases (increases) the "mean" entropy in the regions 
with upward (downward) local flows, thus diminishing 
the buoyancy forces and reducing the "mean" vertical 
velocity Uz and the "mean" vorticity W. This mecha- 
nism dampens the large-scale instability. 

The above two competitive effects determine the 
growth rate of the large-scale instability. Perturbation 
analysis of the mean-field equations in the shear-free con- 
vection regime results in the following expression for the 
growth rate 7 of long- wave perturbations: 



7 (xgF^r'^K^ sin 9\ 



3 5a 2 

a — — — sm 

8 4 



-,1/2 




(4) 



FIG. 2: Experimental set-up: (1) - walls of the chamber; (2) - 
vertically oriented stationary grid; (3) - cooled top with a heat 
exchanger; (4) - laser light sheet (5) - heated bottom with a 
heat exchanger; (6) - two progressive-scan 12 bit digital CCD 
cameras. 



complex but " mean" flow. The growth rate of the large- 
scale instability depends also on the degree of anisotropy 
of turbulent velocity field. The detailed analysis of this 
facet is performed in Section 4 in order to explain the 
obtained experimental results. 



III. EXPERIMENTAL SET-UP 



(Elperin et al. 2002; 2005), where the parameter f3 = 
(I K)~^ 3> 1, I is the maximum scale of turbulent mo- 
tions, 6 is the angle between the vertical unit vector e 
and the wave vector K of small perturbations, g is accel- 
eration of gravity. Equation was derived for the case 
N'^ ^ gF*TK^, where N is the Brunt- Vaisala frequency. 
The large-scale instability occurs when a(5cos^ 9 — 1) > 
3/2. This yields two ranges for the instability: 



2(5cos2 6i - 1; 



< a < 3 



9 3 

< a < 

2 2(1 



5 cos^ 



(5) 
(6) 



where we took into account that the parameter a varies 
in the interval —9/2 < a < 3. The first range for 
the large-scale instability in Eq. ||SJ) is for the angles 
3/10 < cos^ e <1 (the ratio < L^/Lh < 1.53), and the 
second range for the large-scale instability (see Eq. ©) 
corresponds to the angles < cos^ 6* < 2/15 (the ratio 
2.55 < Lz/Lh < 00), where L^/Lh = Kh/K^ = ianO 
and Kh is the horizontal component of the wave vec- 
tor. The maximum growth rate of the large-scale insta- 
bility is achieved at the scale of perturbations L « 10 Z, 

where L = \/ ^ Lz"^ + . The characteristic time of 
excitation of the large-scale instability is of the order of 
20 — 30 T. Therefore, the typical length and time scales 
of these structures are much larger than the turbulence 
scales. This justifies the separation of scales assump- 
tion required for treating coherent structures as albeit 



Now we describe the experimental study of large-scale 
circulations formed in turbulent convection. The exper- 
iments were conducted in a rectangular chamber with 
dimensions 26 x 58 x 26 cm in air flow (see Fig. 2). The 
side walls of the chamber are made of transparent Per- 
spex with the thickness of 10 mm. A number of experi- 
ments were also conducted with different additional ther- 
mal insulation of the side walls of the chamber in order to 
study whether a heat flux through the side walls affects 
the turbulent convective pattern. First, the side walls of 
the chamber were insulated with Styrofoam plates with 
low thermal conductivity (k ^ 0.033 W/ (m K)) and with 
the thickness of 30 mm. Two of the side plates were re- 
moved for a short time when the images of the flow were 
recorded. In the next series of experiments we installed 
additional Perspex plates with a thickness of 6 mm and 
with an air gap of 2 mm between these plates and the 
outside walls of the chamber. Finally, we performed ex- 
periments where these two types of thermal insulation 
were used simultaneously. 

One vertically oriented stationary grid (with 50 mm 
mesh and 10 mm bars where bars are arranged in a square 
array) was attached to the left horizontal rod. This grid 
was positioned at a distance of one mesh size from the 
left wall of the chamber parallel to it. Hereafter, we use 
the following system of coordinates: Z is the vertical axis, 
the F-axis is perpendicular to the grid and the XZ-plane 
is parallel to the grid plane. 

A vertical mean temperature gradient in the turbulent 
flow was formed by attaching two aluminium heat ex- 
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changers to the bottom and top waUs of the test section 
(a heated bottom and a cooled top wall of the chamber). 
The thickness of the aluminium plates is 2.5 cm. The 
top plate is a bottom wall of the tank with cooling wa- 
ter. Temperature of water circulating through the tank 
and the chiller is kept constant within 0.1 K. Cold wa- 
ter is pumped into the tank through two inlets and flows 
out through two outlets located at the side wall of the 
chamber. The bottom plate is attached to the electrical 
heater that provides constant and uniform heating. The 
voltage of a stable power supply applied to the heater 
varies in the range from 35 V to 200 V. The power of the 
heater varies in the range from 10 W to 300 W. The tem- 
peratures of the conducting plates were measured with 
two thermocouples attached at the surface of each plate. 
The temperature difference between the top and bottom 
plates, AT, varies in the range from 5 K to 80 K depend- 
ing on the power of the heater (i.e., the Rayleigh number 
was changed within the range of Ra = 10^ — 1.6 x 10^). 
The temperature in the probed region was measured with 
a specially designed temperature probe with twelve sen- 
sitive thermocouples. The temperature measurements 
showed that the thermal structure inside the large-scale 
circulation is inhomogeneous and anisotropic. The hot 
thermal plumes accumulate at one side of the large-scale 
circulation, and cold plumes concentrate at the opposite 
side of the large-scale circulation. 

The velocity field was measured using a Particle Im- 
age Velocimetry (PIV), sec Raffel ct al. (1998). A digi- 
tal PIV system with LaVision Flow Master III was used. 
A double-pulsed light sheet was provided by a Nd-YAG 
laser (Continuum Surelite 2 x 170 mJ). The light sheet 
optics includes spherical and cylindrical Galilei telescopes 
with tuneable divergence and adjustable focus length. 
We used two progressive-scan 12 bit digital CCD cam- 
eras (with pixel size 6.7 /urn x6.7/Ltm and 1280 x 1024 
pixels) with a dual-frame-techniquc for cross-correlation 
processing of captured images. A programmable Timing 
Unit (PC interface card) generated sequences of pulses 
to control the laser, camera and data acquisition rate. 
The software package LaVision DaVis 7 was applied to 
control all hardware components and for 32 bit image ac- 
quisition and visualization. This software package com- 
prises PIV software for calculating the flow fields using 
cross-correlation analysis. 

An incense smoke with sub-micron particles as a tracer 
was used for the PIV measurements. Smoke was pro- 
duced by high temperature sublimation of solid incense 
particles. Analysis of smoke particles using a microscope 
(Nikon, Epiphot with an amplification of 560) and a PM- 
300 portable laser particulate analyzer showed that these 
particles have an approximately spherical shape and that 
their mean diameter is of the order of 0.7/Ltm. The proba- 
bility density function of the particle size measured with 
the PM300 particulate analyzer was independent of the 
location in the flow for incense particle size of 0.5 — 1 /xm. 

Scries of 130 pairs of images acquired with a frequency 
of 4 Hz were stored for calculating the velocity maps and 



FIG. 3: Mean flow patterns during increase of the the tem- 
perature difference between the bottom and the top walls of 
the chamber: (a) AT = 10 K (Ra = 0.18 x 10*); (b) AT = 20 
K (Ra = 0.33 x 10^); (c) AT = 70 K (Ra = 0.8 x 10*). 



FIG. 4: Mean flow patterns during decrease of the the tem- 
perature difference between the bottom and the top waUs of 
the chamber: (a) AT = 60 K (Ra = 0.74 x 10*); (b) AT = 35 
K (Ra = 0.52 x 10*); (c) AT = 15 K (Ra = 0.26 x 10*). 

for ensemble and spatial averaging of turbulence charac- 
teristics. The center of the probed flow region coincides 
with the center of the chamber. We measured the velocity 
field in the flow areas 487.5 x 240 mm^ with a spatial res- 
olution of 234.7/Ltm / pixel. These regions were analyzed 
with interrogation windows of 32 x 32 pixels. A velocity 
vector was determined in every interrogation window, al- 
lowing us to construct a velocity map comprising 65 x 32 
vectors. 

Mean and r.m.s. velocities, two-point correlation fimc- 
tions and an integral scale of turbulence were deter- 
mined from the measured velocity fields. The mean and 
r.m.s. velocities for each point of the velocity map (1024 
points) were determined by averaging over 130 indepen- 
dent maps. The two-point correlation functions of the 
velocity field were determined for each point of the cen- 
tral part of the velocity map (32 x 32 vectors) by aver- 
aging over 130 independent velocity maps. An integral 
scale I of turbulence was determined from the two-point 
correlation functions of the velocity field. These measure- 
ments were repeated for various temperature differences 
between the bottom and the top walls of the chamber. 
The size of the probed region did not affect our results. 
Similar experimental set-up and data processing proce- 
dure were used in experimental study by Eidelman et al. 
(2002; 2004) and Buchholz et al. (2004) of a new phe- 
nomenon of turbulent thermal diffusion (see Elperin et 
al. 1996; 1997; 2001). 

The maximum tracer particle displacement in the ex- 
periment was of the order of 8 pixels, i.e., 1/4 of the in- 
terrogation window. The average displacement of tracer 
particles was of the order of 2.5 pixels. Therefore, the av- 
erage accuracy of the velocity measurements was of the 
order of 4% for the accuracy of the correlation peak de- 
tection in the interrogation window of the order of 0.1 
pixels (see, e.g., Adrian 1991; Westerweel 1997; 2000). 



IV. HYSTERESIS PHENOMENON 

In this Section we discuss the hysteresis phenomenon 
which was found in the turbulent convection by varying 
the temperature difference between the bottom and the 
top walls of the chamber within a range of 5 to 80 K 
(the Rayleigh number was changed within the range of 
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Ra w 10^ — 10^). The hysteresis loop comprises the two- 
cell and one-cell flow patterns while the aspect ratio of 
the chamber A = 2.23 is kept constant. In particnlar, 
increasing the temperature difference from 5 to 15 K we 
observed first two-cell flow pattern with the downward 
motions in the central region of the chamber between 
two cells (Fig. 3a), then one-cell flow pattern within a 
range of the temperature difference of 20 to 30 K with the 
counterclockwise mean flow (Fig. 3b). Further increase 
of the temperature difference from 35 to 80 K results in 
two-cell flow pattern with the upward mean flow in the 
central region of the chamber between the two cells (Fig. 
3c). Decreasing the temperature difference we observed 
two-cell flow pattern within a range from 80 to 45 K with 
the upward flow in the region between two cells (Fig. 
4a), then one-cell flow pattern with clockwise mean flow 
within a range of the temperature difference of 40 to 20 K 
(Fig. 4b). Further decrease of the temperature difference 
from 15 to 5 K results in two-cell flow pattern with the 
downward mean flow in the central region of the chamber 
between the two cells (Fig. 4c). 

During the transition from the two-cell to one-cell flow 
patterns, the direction of the mean flow in the one-cell 
pattern coincided with that of the right cell in the two-cell 
flow pattern. The latter may be caused by an asymmetry 
of the chamber due to the presence of the grid near the 
left wall. This grid introduces an additional anisotropy in 
the small-scale tiirbulent convection. In all experiments 
the transition from the two-cell to one-cell flow patterns 
was found to be accompanied by the change of the sign 
of the degree of anisotropy x of turbulent velocity field, 
which is defined as 

X=^[^-ll. (7) 

In particular, the dependence of the degree of anisotropy 
X of turbulent velocity field on the Rayleigli munber ob- 
tained in the experiment is shown in Fig. 5. The param- 
eter X was negative for the two-cell flow pattern and was 
positive for the one-cell flow pattern. 

These experimental observations can be explained by 
invoking the theory of formation of coherent structures 
developed by Elperin ct al. (2002; 2005). In particular. 
Fig. 6 shows the theoretical dependence of the growth 
rate of the large-scale instability versus the parameter x 
for the one-cell and two-cell flow patterns. Note that the 
aspect ratio of the large-scale cell in the two-cell flow pat- 
tern observed in the experiment was approximately 1 (see 
Figs. 3a; 3c; 4a; 4c), while in the one-cell flow pattern it 
was approximately 2 (see Figs. 3b; 4b). The growth rate 
of the large-scale instability for the two-cell mode, 72, is 
larger than that for the one-cell mode, 71, for negative 
values of the degree of anisotropy x of turbulent velocity 
field, and 72 < 71 for positive values of the parameter x- 
Figure 6 is plotted for the degree of thermal anisotropy 
a = 1.7. 

The difference A7 = 72—71 of the growth rates of 
the large-scale instability for the two-cell and one-cell 
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FIG. 5: Degree of anisotropy x of turbulent velocity field 
versus the Rayleigh number: increase of Ra (solid line) and 
decrease of Ra (dashed line). 




FIG. 6: Growth rate of the large-scale instability versus the 
parameter x for the one-cell flow pattern (dashed line) and 
for the two-cell flow pattern (solid line), for a = 1.7 and 
L/lo = 10. 



flow patterns versus the degree of thermal anisotropy 
a is shown in Fig. 7. The range of parameters where 
A7 > 0, corresponds to the two-cell flow pattern, since 
in this range the growth rate of the excitation of the two- 
cell mode is larger than that of the one-cell mode. Ac- 
cordingly, the range of parameters where A7 < 0, corre- 
sponds to the one-cell flow pattern, because in this range 
the growth rate of the excitation of the two-cell mode is 
smaller than that of the one-cell mode. Let us deflne the 
point X = Xb for a given parameter a where A7 = 0. The 
dependence of the bifurcation point Xb versus the degree 
of thermal anisotropy a is shown in Fig. 8. The value Xb 
increases with the increase of the parameter a. 

Note that the range of the degree of thermal anisotropy 
a when two kinds of flow patterns (the two-cell and one- 
cell flow patterns) can exist in the system is very nar- 
row. For instance, when —0.3 < X ^ 0.3 the parameter 
a changes within the range between 1.63 and 1.77 (see 
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FIG. 7: The difference A7 = 72 — 71 in the growth rate 
of the large-scale instability for the two-cell and one-coll flow 
patterns versus the degree of thermal anisotropy ol for differ- 
ent values of the parameter x- X~ —0.3 (thick curve) and 
X = 0.3 (thin curve), for L//o = 10. The range with A7 > 
corresponds to the two-cell flow pattern (dashed line), and 
the range with A7 < is for the one-cell flow pattern (solid 
line). 



gree of anisotropy x of turbulent velocity field versus the 
Rayleigh number for the chamber with additional ther- 
mal insolation with Perspex plates (dashed line) and for 
the experiments without the additional thermal insola- 
tion (solid line) in the case of increase of the Rayleigh 
number (Fig. 9a) and during decrease of Ra (Fig. 9b). 
The difference in the dependence x(R'a) was not essential, 
and the hysteresis phenomenon was observed in these ex- 
periments as well. 

Notably, the vertical grid (see Fig. 2) which has been 
placed asymmetrically, introduces additional anisotropy 
in turbulent fluid flow in a chamber. Without this ver- 
tical grid the dynamics of flow patterns becomes more 
complicated. In particular, after a two-cell flow pat- 
tern at higher values of the Rayleigh number we often 
observed an additional transition to the one-cell flow 
pattern. However, we did observe the hysteresis phe- 
nomenon in this case. When we also changed the aspect 
ratio from A = 2 to A = 2.23, the hysteresis phenomenon 
was observed as well. 



CONCLUSIONS 
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FIG. 8: Bifurcation point Xb (where A7 = 0) versus the 
degree of thermal anisotropy a for L/io = 10. 



Fig. 7). The latter theoretical prediction is in an agree- 
ment with the experimental observations. In particular, 
our preliminary measurements show that the degree of 
thermal anisotropy a in the chamber in the experiments 
was of the order of 1 — 2. Therefore, the above discus- 
sion shows that the obtained experimental results are in 
a good agreement with the predictions of the theory of 
formation of coherent structures in turbulent convection 
(see Elperin et al. 2002; 2005). 

A number of experiments were also conducted with 
different additional thermal insulation of the side walls 
of the chamber in order to study whether a heat flux 
through the side walls affects the turbulent convective 
pattern. In particular, we plotted in Fig. 9 the de- 



In the present study we demonstrated that the 
anisotropy of turbulent velocity field plays a crucial role 
in the hysteresis phenomenon that was found in our ex- 
periments in turbulent convection. The hysteresis loop 
comprises the two-cell and one-cell flow patterns while 
the aspect ratio of the chamber is kept constant. The 
observed transition from the two-cell to one-cell flow pat- 
terns causes a drastic change of the degree of anisotropy 
X of the turbulent velocity field from negative to pos- 
itive values. This finding is in a good agreement with 
the theoretical predictions. In particular, we found that 
there is a critical value Xh of the degree of anisotropy of 
turbulent velocity field when the growth rate of the large- 
scale instability for the one-cell flow pattern is larger than 
that for the two-cell flow pattern. This critical value de- 
pends on the degree of thermal anisotropy. Transition 
through this critical value Xfc results in a bifurcation, 
i.e., it causes the transition from the two-cell to one-cell 
flow patterns. Then further increase of the temperature 
difference between the bottom and the top walls of the 
chamber decreases the degree of anisotropy Xi and the 
system passes again through the bifurcation from one-cell 
to two-cell flow patterns. Decreasing the temperature dif- 
ference causes these bifurcations but for different values 
of the temperature difference. Therefore, we observed 
the hysteresis phenomenon in turbulent convection. The 
developed theory of coherent large-scale circulations in 
turbulent convection (Elperin et al. 2002; 2005) is in a 
good agreement with the experimental observations. 

Note that the nonlinear dynamics may play a very es- 
sential role in the observed transitions between one-cell 
flow pattern to the two-cell flow pattern (and vice versa) 
when we change the values of the Rayleigh numbers. 
After parametrization of the turbulence the nonlinear 
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mean-field cqiiations determine the dynamics of modes 
(flow patterns). In our experiments witli A = 2 — 2.23 we 
have a two-mode system (one-cell and two-cell modes), 
and nonlinear interactions between the modes play a very 
essential role. Experiments and numerical simulations 
with large aspect ratios (Willis et al. 1972; Fitzjarrald 
1976; Hartlep et al. 2003) can be regarded as a multi- 
mode system, which may be strongly different from the 
two-mode system. 

The obtained results might be important in atmo- 
spheric turbulent convection (see, e.g., Zilitinkevich 1991; 
Etling and Brown 1993; Atkinson and Wu Zhang 1996; 
Briimmer 1999) and industrial turbulent flows (see, e.g., 
Kakac et al. 1987; Incropera and DeWitt 2002). 
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FIG. 9: Degree of anisotropy x oi turbulent velocity field 
versus the Rayleigh number for the experiments shown in Fig. 
5 (solid line) and for the chamber with additional thermal 
insolation (dashed line) during: (a) increase of Ra and (b) 
decrease of Ra. 
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